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Abstract 

In the context of phase-space quantization, matrix elements and observables result from integration 
of c-number functions over phase space, with Wigner functions serving as the quasi-probability measure. 
The complete sets of Wigner functions necessary to expand all phase-space functions include off-diagonal 
Wigner functions, which may appear technically involved. Nevertheless, it is shown here that suitable 
generating functions of these complete sets can often be constructed, which are relatively simple, and 
lead to compact evaluations of matrix elements. New features of such generating functions are detailed 
and explored for integer-indexed sets, such as for the harmonic oscillator, as well as continuously indexed 
ones, such as for the linear potential and the Liouville potential. The utility of such generating functions 
is illustrated in the computation of star functions, spectra, and perturbation theory in phase space. 



1 Introduction 



General phase-space functions f(x,p) and g(x,p) compose noncommutatively through Groenewold's *- 
product |l|], which is the unique associative pseudodifferential deformation |Q| of ordinary products: 

^_ e ih(d x dp- dpdx)/2 _ /-Q 

This product is the cornerstone of deformation (phase-space) quantization ||, |||, |], || , as well as applications 
of matrix models and non-commutative geometry ideas in M-physics Q. Its mechanics, however, is not 
always straightforward. 

The practical Fourier representation of this product as an integral kernel has been utilized widely since 
Baker's early work, 

/ * g = _L J dp' dp" dx' dx" f{x',p') g(x",p") exp (p(x' - x") + p'(x" - x) + p"(x - x'))^j . (2) 
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The determinantal nature of the star product controls the properties of the phase-space trace || [|, 

J dpdx f* 9 = J dpdx f 9 = J dpdxg *f. (3) 

The above ★-product and phase-space integrals provide the multiplication law and, respectively, the 
trace in phase-space quantization the third autonomous and logically complete formulation of quan- 
tum mechanics beyond the conventional formulations based on operators in Hilbert space or path inte- 
grals. (This formulation is reviewed in J2], ||.) Properly ordered operators (e.g., Weyl-ordered) correspond 
uniquely to phase-space c-number functions (referred to as "classical kernels" of the operators in ques- 
tion); operator products correspond to ★-products of their classical kernels; and operator matrix elements, 
conventionally consisting of traces thereof with the density matrix, correspond to phase-space integrals of 
the classical kernels with the Wigner function (WF), the Weyl correspondent of the density matrix |IQ| , ||]. 
The celebrated ★-genvalue functional equations determining the Wigner functions [^, [I!]] and their spectral 
properties (e.g. projective orthogonality |]l2|) are reviewed and illustrated in ||. 



The functions introduced by Wigner [10] and Szilard correspond to diagonal elements of the density 
matrix, but quantum mechanical applications (such as perturbation theory), as well as applications in 
noncommutative soliton problems [13 often require the evaluation of off-diagonal matrix elements; they 
therefore utilize the complete set of diagonal and off-diagonal generalized Wigner functions introduced 
by Moyal ||. For instance, in noncommutative soliton theory, the diagonal WFs are only complete for 
radial phase-space functions (functions ★-commuting with the harmonic oscillator hamiltonian — the radius 
squared), whereas deviations from radial symmetry necessitate the complete off diagonal set. 



As for any representation problem, the particular features of the ★-equations under consideration fre- 
quently favor an optimal basis of WFs; but, even in the case of the oscillator, the equations are technically 
demanding. It is pointed out here, however, that suitable generating functions for them, acting as a trans- 
form of these basis sets, often result in substantially simpler and more compact objects, which are much 
easier to use, manipulate, and intuit. Below, after some elementary overview of the Weyl correspondence 
formalism (Sec 2), we illustrate such functions for the harmonic oscillator (Sec 3), which serves as the 
archetype of WF bases indexed discretely; it turns out that these generating functions amount to the 
phase-space coherent states for WFs, and also the WFs of coherent state wavefunctions (Appendix A). 
Direct applications to first order perturbation theory are illustrated in Appendix B. 

For sets indexed continuously, the generating function may range from a mere Fourier transform, 
illustrated by the linear potential (Sec 4), to a less trivial continuous transform we provide for the Liouville 
potential problem (Sec 5), where the advantage of the transform method comes to cogent evidence. 

Throughout our discussion, we provide the typical ★-composition laws of such generating functions, 
as well as applications such as the evaluation of ★-exponentials of phase-space functions (Appendix C), 
or ★-versions of modified Bessel functions (techical aspects of integral transforms of which are detailed in 
Appendix D). Appendix E provides the operator (Weyl-) correspondent to the generating function for the 
Liouville diagonal WF introduced in Sec 5. 
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2 Overview of General Relations in the Weyl representation 

Without loss of generality, we review basic concepts in two-dimensional phase space, (x,p), as the extension 
to higher dimensions is straightforward. In addition, we first address discrete spectra, E n , n = 0, 1, 2, 3, 
and will only later generalize to continuous spectra. 

In the Weyl correspondence [14|, c-number phase-space kernels a(x,p) of suitably ordered operators 
A(X ,V) are defined by 

a(*,p)^ J dye-^(x-^y\A(X,V)\x + ^y). (4) 
Conversely, the ordering of these operators is specified through 

A(X,V) = 7 — / drdadxdp a(x, p) exp(ir(V — p) + ia(X — x)). (5) 
(2tiT J 

An operator product then corresponds to a star-composition of these kernels [|l]] , 

1 f ft ft 

a(x,p)*b(x,p) = — I dye- iyp (x--y\A(X,V)B(X,V)\x + -y), (6) 

Moyal H appreciated that the density matrix in this phase-space representation is a hermitean gener- 
alization of the Wigner function: 



f mn (x,p) = J dy e iyp {x - ^y\tp n ) {ipm\x + ^y) 



= ^~ / ' dye WP ^*m{ X ~ \y)^n{x + ~y) = f* m (x,p) , (7) 

where the tp m (x)s are (ortho-)normalized solutions of a Schrodinger problem. (Wigner [10] mainly con- 
sidered the diagonal elements of the density matrix (pure states), usually denoted as f m = / mm .) As a 
consequence, matrix elements of operators are produced by mere phase space integrals ||, 

{ipm\A\ip n ) = J dxdp a(x,p)f mn (x,p). (8) 

The standard machinery of density matrices then is readily transcribed in this language, e.g. the trace 
relation ||, 

J dxdp fmn(x,p) = J dx t/j*(x)t/j m (x) = 5 m n | (9) 

and § 

fmn * fkl = "Z T^mlfkn = "T^rnlfkn ■ (10) 

2irn h 

Given (||), it follows from the above that Q 

dxdp f mn (x,p)f? k (x,p) = — <5 m ^ n jfe ■ (H) 

For complete sets of input wavefunctions, it also follows that ||: 

^2fmn(x,P)fmn( x ',P') = ^$(x ~ x')5(p - p') . (12) 
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An arbitrary phase-space function tp(x,p) can thus be expanded as 

tp(x,p) = ^c mn / m „(x,p), (13) 



the coefficients being specified through (|11|), 

c mn = 2irh I dxdp f^ n (x,p)(f(x,p). (14) 



Further note the resolution of the identity 0] , 



Ef^p) = ^n = l- ( 15 ) 



For instance, for eigenfunctions of the hamiltonian 7i(X,V) with eigenvalues E n , the corresponding 
WFs satisfy the following star-genvalue equations || (further cf (ll], with H(x,p), the phase-space 
kernel of H(X,V): 

H * fmn — E n f mn , fmn * H — E m f mn- (16) 

The time dependence of a pure state WF is given by Moyal's dynamical equation ||: 

ih^- t f(x,p;t)=H*f{x,p;t)-f(x,p;t)*H. (17) 
By virtue of the *- unitary evolution operator (a "^-exponential" ||]), 

U+ix, p; t) = el tH/h = 1 + (it/h)H(x, p) + * H + * H * H + ... , (18) 

the time-evolved WF is obtained formally in terms of the WF at t = 0, 

f(x,p;t) = U- 1 (x,p;t)*f(x,p;0)*U*{x,p;t). (19) 

(These associative combinatoric operations completely parallel those of operators in the conventional for- 
mulation of quantum mechanics in Hilbert space [15|.) Just like any star- function of H, this ★-exponential 



can be computed, 16] 



exp^itH/h) = exp*(itH/h) ★ 1 = exp^itH/h) * 2vrft^ f nn = 2irh^2e ltEn/h f nn . (20) 

n n 

(Of course, for t = 0, the obvious identity resolution is recovered). 

For continuous spectra, the sums in the above relations extend to integrals over a continuous parameter 
(the energy), and the Kronecker 5 mn s into 5-functions (these last ones reflecting the infinite normalizations 
of unnormalizable states). E.g., eqns (||, |ll|) extend to 



J dxdp f El E 2 (x,p) = 5(E 1 - E 2 ), (21) 
J dxdp f El E 2 {x,p)r E , E ,{x,p) = 2^(#i - E[)5(E 2 - E' 2 ), (22) 
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Completeness (|12|) extends to 

J dE l dE 2 f El E 2 (x,p)r ElE2 (x',p') = J^5(x-x')5{p-p'). (23) 
More generally, ([To|) extends to 

f El E 2 * f E >E> = ^>(Ei ~ E' 2 ) f E[E2 . (24) 

Finally, eqn (|i~5l) extends to 

1 1 f J.. _-imi f . , % 



2tt7i 2vr 

d£ f EE {x,p), (25) 



and hence (^) extends to 

exp* (itH/h) = 2nh J dE e itE/h f EE (x, p) . (26) 

3 Generating Functions for the Harmonic Oscillator 

Consider the harmonic oscillator, 

H(x,p) = l(p 2 + x 2 ), (27) 

where, without loss of generality, parameters have been absorbed in the phase space variables: m = 1, 
ui = 1. Further recall that the normalized eigenfunctions of the corresponding operator hamiltonian TC are 

_ 1 12 

ip n (x) = (\/7r2 n n\) 2 e~z x H n (x), for the eigenvalues E n = h(n + 1/2). Define a radial and an angular 
variable, 

z = 4H = 2(x 2 + p 2 ) , tan0 = ~, (28) 

x 

so that 

i 

aV2 =(x + ip) = \x + ip\e ld = 2 e ie . (29) 



Groenewold |T|, as well as Bartlett and Moyal [17], have worked out the complete sets of solutions 



to Moyal's time-evolution equation (|17|), which are all linear combinations of terms exp(ii(m — n)) f mn . 
They solved that equation indirectly, by evaluating the integrals (^) for time-dependent Hermite wavefunc- 
tions, which yield generalized Laguerre polynomial-based functions. More directly, Fairlie || dramatically 
simplified the derivation of the solution by relying on his fundamental equation (|l6|). He thus confirmed 
Groenewold's WFs fill, Il7|, 



f mn (x, P ) = t})lM z ^ e -^(n- m) e L n m{ ^ m) 

it V n\ 



The special case of diagonal elements, 



(—l) n 

fn = fnn= K -^e- Z l 2 L n {z), (31) 
7T 
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constitutes the time-independent "★-genfunctions" of the oscillator hamiltonian kernel j| (i.e. the complete 
set of solutions of the time-independent Moyal equation H-kf—f-kH = 0, where H-kf n = E n f n . Incidentally, 
(|l0| ) restricted to diagonal WFs closes them under ★-multiplication fli~2|] , f m ★ /„ = 5 mn f m /(27rh).) That 
is to say, "radially symmetric" phase-space functions, i.e. functions that only depend on z but not 9, 
can be expanded in terms of merely these diagonal elements — unlike the most general functions in phase 
space which require the entire set of off-diagonal f mn above for a complete basis. Note, however, that all 
★-products of such radially symmetric functions are commutative, since, manifestly, 

^ Cnfn * ^ d m f m — ^ d m f m ★ ^ Cnfn- (32) 



Moreover, the ★-exponential (^) for this set of ★-genfunctions is directly seen to amount to 

exp*(itH/h) = ^cos(^ exp ^Ftan(^ , (33) 

which is, to say, a Gaussian in phase space 0. As an application, note that the hyperbolic tangent ★- 
composition law of gaussians follows trivially (since these amount to ★-exponentials with additive time 
intervals, exp*(t/) ★ exp Jr (T/) = exp^((t + T)f), j§, 

1 / a + b , 2 2 , 



exp(--(z^)j ★ exp^-^+^j =_ ex ^- w —- ) ^ +p ^ ] . (34) 



We now introduce the following generating function for the entire set of generalized Wigner functions, 



G(a,(3;x,p) = V _ ^ fmn 



-z-ie- w a) L n ~ m {z) 



(35) 



Utilizing the identity pf 8.975.2, 



]T L^ m {z) k m = e- zk (l + k) n , 



m=0 



we obtain 



z 2 e l9 a 



7r n\ \ 



ze ie -a(3) e 



Thus, 



7T 



1 / • Z 

G(a, /?; x,p) = — exp I \/z(«e~ ! + fie 1 ) — a/3 — - 
7T V 2 



(36) 



(37) 



(38) 
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Since 



y/z(ae~ ie + pe ie ) = V2(a + (3)x - s/2ip(a - 0), 



one can re-express: 



G(a, (3; x,p) = G*(/3, a; x,p) = — exp (a/3 — ( x — — 

tt \ \ V2 

As the name implies, from G(a,/3;x,p), the / mn s are generated by 

1 d m d n 



p + i 



a — (3 

~7T 



fmn{%i P) 



G(a,P;x,p) 



a=/3=0 



(39) 



(40) 



(41) 



These functions ^-compose as 



The phase-space trace is 



G(a,[3)*G(e,() 



2Tih 



j dxdp G(a, (3) 



G(e,(3). 



By ( pi ) , the action of the Hamiltonian kernel on this function is 

H*G = h(± + (3^j G = h( K ~-a{3 + /3^e i0 ) G, 



and 



Consequently, 



G * II - h ( i + a-^J G = h(^-a(3 + aj~ze~ ie ) (V . 



J dxdp H * G(a, /3) = h Q + /3 ^ e a/3 = Q + a/3 J e 

The spectrum then follows by operating on both sides of this equation, 

1 d n d n 



ai3 



En 



nl da n d(3 
h d n d n (I 



dxdp H ★ G(a, (3) 
+ a/3) e af3 



a=/3=0 



a=/3=0 



h\- + n 



(42) 
(43) 

(44) 
(45) 
(46) 



In general, matrix elements of operators may be summarized compactly through this generating 
in phase space. 



(47) 
function 



This generating function could be interpreted as a phase-space coherent state, or the off-diagonal WF 
of coherent states, as discussed in Appendix A, [HI 



G(a,(3;x,p) = exp„(/3a T ) f exp„(aa) 



(48) 
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a*G(a,p) = H0G(a,p), a f * G(a, p) = -^G(a, p) 

op 

d 

G(a,P)*a = T^G(a,P), G(a, 0) * a f = haG(a, p), (49) 
and hence eqs (|S^j45| ) amount to 

H*G(a,P) = (aUa + ^)*G(a,P)=h(p-^ + ^\G( a ,P) 
G(a,P)*H = G(a,p)*(ai*a+£) = H(a^ + ^\G(a,p). (50) 
This formalism finds application in, e.g., perturbation theory in phase space, cf. Appendix B. 

4 Generating Functions for the Linear Potential 

The linear potential in phase space has been addressed [|ll|] (also see [19, We shall adopt the simplified 



conventions of Q, ie. m = 1/2, h = 1. The Hamiltonian kernel is then, 

H(x,p) = p 2 + x , (51) 
and the eigenfunctions of 7i are Airy functions, 



i r+oo 

$ E {x) = — dX e iX(z-*-x 2 /3) = Ai ( x _ ; ( 52 ) 

indexed by the continuous energy E. The spectrum being continuous, the Airy functions are not square 
integrable, but have continuum normalization, f dxifrg (x)ipE 2 ( x ) = $ {E\ — E2), instead. Thus, (21) et 
seq. are now operative. The generalized WFs are p} 



= e ^- E ^M^/ 3 (x+p 2 -^±^)y (53) 
The ★-exponential ( |26| ) then is again a plain exponential of the shifted hamiltonian kernel, 

exp, (it (x + p 2 ))= 2vr y°° dEe iEt ^ Ai (V/ 3 (x + p 2 - ^^ X) = exp [it (x + p 2 + t 2 /U)) . 

(54) 

(This could also be derived directly, as the CBH expansion simplifies dramatically in this case, cf. Appendix 
C.) As before, the ^-composition law for plain exponentials of the hamiltonian kernel function follows, 

exp (a (x + p 2 )) * exp (6 (x + p 2 )) = exp ^(a + b) ^x + p 2 — ~^ a ^j ^ • (55) 



Since the complete basis Wigner functions are now indexed continuously, a generating function from 
them must rely on an integral instead of an infinite sum. The simplest transform is possibly a double Fourier 
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transform with respect to the energy indices (but note the transform factors exp(iE\X), exp(—iE 2 Y) may 
also be regarded as plane waves). Suitably normalized, 



f+oo r+oo 

G(X,Y;x,p) = 2tt / dE 1 \ dE 2 



JElX 



2vr 



Ie 1 e 2 (x,p) 



1 



-iE 2 Y 



2vr 



(56) 



-L | + °° dE X J +CO dE 2 e i(E 1 -E 2 ) P+ iE 1 X-iE 2 Y 2 2/3 Ai ^2/3 ^ +p 2 _ ^^X\ 

oo . 

doj e ^P+i(E+u/2)X-i(E-u/2)Y 2 2/3 Ai ^ + p 2 _ ^ 

o 

d£ e ^(^- y ) 2 2 / 3 Ai (2 2 / 3 (x +p 2 - B) 

J —oo 



i r+oo f+oo 

— / d£ 

2vr 



5lp + 



X + Y 
2 

x + y 

2 

X + Y 



oo 

+oo 1 

dEe iE(X-Y) J_ / 

-oo 2?r 

,i(x-y)(2:+p 2 +(x-y) 2 /i2) 



The phase-space trace is 



J dxdp G(X, Y; x,p) = 2tt5{X - Y) , 



(57) 



and, given (24) for these functions, /exE 2 * Ie'e' = ^k^^ 1 ~ ^2) ^ ne ^-composition law for these 
Gs is 

G(X,Y;x,p)*G(W,Z;x,p) = 5(X - Z) G{W,Y;x,p). (58) 



5 Generating Functions for the Liouville Potential 

A less trivial system with a continuous spectrum is the Hamiltonian with the Liouville potential, 21]. 
In the conventions of 0], (H = 1, m = 1/2), the Hamiltonian kernel is 

H = p 2 + e 2x , (59) 

and the eigenfunctions of the corresponding TC are 

iPe(x) = r E (*) = ^smh(TrVE) K m {e x ), (60) 

with continuum normalizations J dxip E (x) ipE 2 ( x ) = $ {E\ — E 2 ). The modified Bessel function (Cf. |22f| , 
Ch VI, §6.22) can be written in the Heine-Schlafli form, 

1 r°° 

K ip (e x ) = - dX exp (-e x coshX + iXp) = K_ ip (e x ). (61) 

2 J-oo 

The non-diagonal WF is then 

f El E 2 (.x,p) = ^J dye- 2 ^^sinh(^/^)A^( e -^)^sinh(vr V / ^)^ v/ ^(e^). (62) 
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This Wigner function amounts to Meijer's G function, 

Ax 



1 

8^3 



sinh(vr y/Ei) sinh(vr \JE 2 ) Gf 4 



16 



ip + i\[E\ ip — iy/Ei —ip + i\fE~2 —ip — i\/E% 



(63) 



Alternatively, the WF may be written as a double integral representation, 
fE(k)E(q)(x,p) = (64) 
= ^ ^sinh (ny/E(k?) sinh J^fW^pj J dXdY e tkX e iqY (^^) "^Mp (eV4 cosh X cosh F) , 

where J5(jfe) = k 2 ,E(q) = q 2 . This is an inverse integral transform, as in the preceding section, of a 
generating function 



G(JT,y ; x,p) = 



dk 



sinh I iryE{k) 



dq 



sinh I ny/E{q) 



e -ikx e - iqY (a . jp) 



2 /coshF 



) #2ip (eV4 cosh X cosh FJ = G* (y, X; x,p) . 



7r \coshX / 

The form and construction of this G are consequences of (BT1), as detailed in Appendix D. 



(65) 



However, the ^-composition law of this particular generating function is not so straightforward. It is 
singular, as a consequence of the general relation (p4|) and the behavior of the integrand in (65) as k, q — > OQ. 



By some contrast to the above, an alternate generating function for just the diagonal WFs, /ee = fE, 
could be defined through the spectral resolution of the * — K function, 



poo 

g(z;x,p) = .^^(e") = 2tt J dE K iyfE (e z ) f E (x,p) . 



This can be evaluated by reliance on Macdonald's trilinear identity 23], 



o 



dE K iVE^ Z )^E(x)r E (y) = l exp ' 



- (e x+y ~ z + e x ~ y+z + e~ x+y+z ^ 



1 The singularity may be controlled by regulating the ^-product through imaginary shifts in the momenta 

1 



C[X,Y;x,p-^\ *G\W,Z;x,p+^ 



, , , - , , i , i cosh Y cosh W , , , 

G(W,Y\x,p)Y(e) e x \ — — ■ (cosh X + cosh Z) 

2tt v ' ' ,1J w \ V cosh X cosh Z v ; 



It follows that one derivative with respect to either of X or Z suffices to eliminate the divergence at e = 0, 

lim d x G (x,Y;x,p- — \ * G (w, Z;x,p + ^ 



1 ^,/,,^,, \ / „ , , T /coshYcoshW. , „ , „, 1 >/l , , r sinhX 

-G (W, Y; x, p) (-ax) In — - — - (cosh X + cosh Z)\= —G (W, Y; x,p) I - tanh X - — - - — - 



(66) 



(67) 



Unlike the situation in (p8|), here the RHS vanishes at X = Z. More symmetrically, 



limdxG^X,r;a;,p-^) * d w G ( W, Z; x,p + | 



sinh X 



2^- G W Y '< *M\ tanhX " cosh X + cosh Z 
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Q then is obtained by replacing x — > x + Y and y — > x — Y, and Fourier transforming by - f dY e 2tpY 
J™dEK iV z(e z )f E (x,p) = ±- J dYe-^ Y exp (~ (e 2 *"* + e*+ 2Y + e*~ 2Y ) 



(68) 



Finally, simplifying the RHS gives 



2n J™dEK iV z(e z ) f E (x,p) = exp(-ie 2 ^ J dY e' 2 ^ exp (~ e z (e 2Y + e~ 2Y ) 



exp ( ~\e 2x ~ z ) A' (> ((■■-) - t/(::.r./0. (09) 



As a side check of this expression, (p9|), note that it must satisfy the equations 

H*g(z;x,p) =G(z;x,p)*H = (-d 2 + e 2z ) G{z; x,p) , (70) 

which follows from the spectral resolution evident in (|66|). Indeed, since e~ z d z Ki p (e z ) = ipe~ z K{ p (e z ) — 
Kip+i (e 2 ), and (— d 2 + e 2z ) Ki p (e z ) = p 2 Ki p (e z ), these relations are satisfied, 

(p 2 + e 2 *)*(exp(-ie 2 — W (**)) = (exp (-1 e 2 ^ ) A* (e 2 )) * (p 2 + e 2x ) 

= exp (-1 e 2x ~ z ^j {-e 2x ' z d z K ip (e z )) + L 2 + X -e 2x ~ z - -e 4 *" 2 ^ exp (~ e 2x ~ z \ K ip (e z ) 

= {-d 2 z + e 2z ) (exp (-1 e 2x ~ z ^ K rp (e 2 )) . (71) 



Parenthetically, as an alternative to the ordinary product form in (p9|), the phase-space kernel Q may 
also be represented as an integral either of a ^-exponential or of a single *-productQ, 

= i J dy exp f~e y ~ z e 2x J * exp(iyp - e z cosh y) . (72) 
This follows from the identities (cf. Appendix C), 



6XP * ' ~ 2sinhy ^ + iyP ) = 6XP V 2^ J * ex P^) = ex P [^2 * + lW ' ' ( '''^ 



The ordinary product form in ( |69| ) and the ★-exponential form in ( [72]) reveal that 0(z; x,p) = G(z; x, —p), 
so one may replace exp (iyp) by cos (yp) in the second line of ( |72[ ) above. Given these, there are several 
ways to verify (|70|). These relations and the star-product expressions for the kernel in ( f72|) are isomorphic 
to those of the corresponding operators, as discussed in Appendix E. 



(74) 



The ^-composition law of these generating functions follows from ( |24[) and Macdonald's identity, 
g(u;x,p)*g(v;x, P ) = ^ J dw exp (~ (e u+v ~ w + e u ~ v+w + e - u+v+w )^j G(w;x,p). 



NB Do not shift the integration parameter y by the phase-space variable x before the star products are evaluated. 
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This also follows directly from the explicit form (|69|) . Again, this is isomorphic to the corresponding 
operator composition law given in Appendix E. 

From the orthogonality of the i^es, the diagonal WFs may be recovered by inverse transformation, 

f smh(iry/~E) 

Je{x,P) = dz — 3 K^{e ) g{z;x,p). (75) 

This representation and the specific factorized x,p-dependence of Q can be of considerable use, e.g., in 
systematically computing diagonal matrix elements in phase space. 

In illustration of the general pattern, consider the first-order energy shift effected by a perturbation 
Hamiltonian kernel H\. It is, cf. Appendix B, eqn fllOj ), 

AE = J dzdxdp E x Smh 2 ^ 3 ) ^( eZ ) 9{zix,p). (76) 

Choosing 
yields 

AE = sinlWg J dzK ^ {eZ) (| dxe2 - exp (_l e ^-^ (JdpK^eF)**^. (78) 



Now, 



dxe 2nx exp ( ~^e 2x - 



2 n - x r (n) e nz , 



and hence, (@ 6.576.4, a = b), 



J dzK^{e z )K vp 



e~ e" 



2 n 3 / n + iy/E + ip \ I n + iy/E — ip \ ( n — iy/E + ip \ ( n — iy/E — ip 
T(n) { 2 J 1 2 J 1 2 J 1 2 



(79) 



(80) 



Thus, 



A£ = ^ — - 4" 2 x 

2vr 3 



( n + — \ / n — iy/E + \ / n — i\AE — 



Finally, (|| 6.422.19), 



4tt G 22 ' ~ s 



T 22 



2-n + i\fE 2-n- iy^ 
2 ' 2 
n J ci\f~E n— i\fE 
2 ' 2 



(81) 



dpe isp/2 F ( n + i ^ +i p \ r [ n + iVE-ip \ ^ ( n - iy/E + ip \ p / n - iy/E - ip \ 



Curtright, Uematsu, & Zachos 



Generating Functions for All WFs 



13 



To sum up, the perturbed energy shift is a Meijer function 

4" sinh(vr/E) 



AE 



8vr 2 



G 22 



22 



2-n + i\fE 2-n- iy/E 
2 ' 2 
n+jy/E n— i\fE 
2 ' 2 



(83) 



In principle, any polynomial perturbation in either x or p can be obtained from this, by differentiation 
with respect to n and s. (Retaining a bit of exponential in x would be helpful to suppress the region of 
large negative x). 
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Appendix A *-Fock Space and Coherent States 

Dirac's Hamiltonian factorization method for algebraic solution of the harmonic oscillator carries 
through (cf. ||) intact in ★-space. Indeed, 

H =~( x -ip)*(x + ip) + ^, (84) 



motivating definition of 

Thus, noting 
and also that, by above, 



a= -y=(x + ip), a) = -j=(x - ip). (85) 

a -k a) — a) * a = H, (86) 

a*f = -Ux + ip)*e-( a ' 2+p2 ) =0, (87) 
v2 

provides a *-Fock vacuum, it is evident that associativity of the ^-product permits the entire ladder 
spectrum generation to go through as usual. The *-genstates of the Hamiltonian, s.t. H * / = / ★ H, are 
thus 

fnn = fn = ^*) n f (*a) n . 



These states are real, like the Gaussian ground state, and are thus left-right symmetric *-genstates. They 
are also transparently ^-orthogonal for different eigenvalues; and they project to themselves, as they should, 
since the Gaussian ground state does, /o * /o = /o/27r^- 

The complete set of generalized WFs can thus be written as 

/ mn = -^=(a t *) n /o(*o) m > rn,n = 0,1, 2, 3, • • • (89) 

yn! ml 
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The standard combinatoric features of conventional Fock space apply separately to left and right (its 
adjoint) ^-multiplication: 

a ★ f n = a* f nn = %y/nf n ,n-\ 

a* * a * f n = h^/na^ k / n>n _i = hnf n 

a) -k f n = a) -k f nn = Vn+ l/ n , n +l 

a -k a) * f n = Vn + l a * f n ,n+i = K n + l)/n, (90) 



f n *a= Vn+ lf n +l,n 

f n *a*a* = h(n + l)/ n 

/n *« f = hy/nf n -l,n 

f n -ka) ka = hnfn ■ (91) 



Furthermore, a left/right (non-self-adjoint) coherent state is naturally defined [|l9|, [24] 
$(a, /3) = exp^(aa t ) / exp^(/5a), a ★ $(a, /3) = a$(a, /3), /3) * a* = /3<I>(a, /?). (92) 



Up to a factor of exp((|a| 2 + |/3| 2 )/2), this is also the WF of coherent states \a) and (/3\, p4| . As indicated 
in the text, this coherent state is identifiable with the generating function G for the harmonic oscillator. 

Appendix B Stationary Perturbation Theory 

Perturbation theory could be carried out in Hilbert space and its resulting wavefunctions utilized to 
evaluate the corresponding WF integrals. However, in the spirit of logical autonomy of Moyal's formulation 
of Quantum Mechanics in phase space, the perturbed Wigner functions may also be computed ab initio 



in phase space [17, 25], without reference to the conventional Hilbert space formulation. The basics are 



summarized below. 

As usual, the Hamiltonian kernel decomposes into free and perturbed parts, 

H = H + \H 1 . (93) 
Fairlie's stationary, real, *-genvalue equations [||, |j for the full hamiltonian, 

H(x,p) kf n {x,p) = f n (x,p)*H(x,p) = E n (\)f n (x,p), (94) 
are solved upon expansion of their components E and / in powers of A, the perturbation strength, 

E n = E n + \En + \ 2 E 2 n + --- (95) 

fn= /„° + A/i + A 2 / 2 + ---. (96) 

Note the superscripts on E and / are order indices and not exponents. Resolution into individual powers 
of A yields the real equations: 

Ho* fn = fn* H = Enfn (97) 
Hq * fn + Hi * fn = fn k H + fn * Hi = E n fn + E n fn (98) 
Hq* fn + Hi* fn = fn* H + f^* Hi = E n fn + E n fn + E^In ■ (99) 
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Left multiplication of (pgj) by yields 

/° * i? * f n + /n * Hi * /° = * + -^n/n * /!. ( 10 °) 

and, by ©, 

f n *H x *f n = Elf n *f n ; (101) 
by ([ll], |l^), and the cyclicity of the trace @, 

J dxdp E\f n *fn = J dxdp(f n *H 1 *f°) = J dxdpiH, * /° * /°) = JL y dsdpHi * /° . (102) 



Hence, 



£?i = y dxdp ffi/o, (103) 



the diagonal element of the perturbation. For the off-diagonal elements, similarly left-*-multiply ( |98| ) by 

J mi 

fm* H 0* fn + fm* H l* fn = E nfm * fn + E nfm * fn • ( 104 ) 

By completeness, f n , i ^ 0, resolves to 
the reality condition dictating 



^n.fcZ ~~ u n,lk 

Consequently, by 10, 



(106) 



E m E a n,klfm * fkl + fm * ^1 * fn ~ E n E a li,klfm * fkl + E ^2^h^ mn ' (^) 
k,l k,l 

and hence 

E m E a n,kmfkm + ^^fm * H\ * /° = ^ a n,kmfkm + E nfn^mn ■ (108) 

For m / n, 

i 

so that 

w o _ 27Tft(/° 



n m 
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Finally, use of (|ll[) , yields 



I. _ ,.i-f,2 / ; . j„ ./','/'./ * fm * II \ * ./', 



(2vr/i) 2 y 



P 2<Kh E° n -E* m 



(2wh) I dxdp Hl * f °* f ™ 1 



E° - E° 



dxdp B x f mn (m^n). (Ill) 



We also have the similar equation for I ^ n. Consequently, a* im is proportional to the matrix element of 
the perturbation, and it vanishes unless I or m is equal to n. (NB This differs from [25] eqn (45).) To sum 
up, 

fn = E f/nm( / dx' dp' H X (x',p')f m n (x',p')) + /°„( / dx>' ^ (x', pO/nrnO*', P'))) ■ (H2) 

By (||), it can be seen that the same result may also follow from evaluation of the WF integrals of perturbed 
wavefunctions obtained in standard perturbation theory in Hilbert space. 

For example, consider Hi = y/2 x = a + a) . It follows that Eq = 0, and 

«n,Zm = Fg0~g0 ) J J dxd Pfmn * ( a + flt ) 



6. 



n,i 



J J dxdp(Vm + 1 /° +1>n + v^TT 



= $n,l (Vm + 1 S m +i, n - Vn + 1 5 mtn+ i ) , (113) 
for m ^ n and the (m «-> I) expression for I ^ n. Hence, 

fn = Mfn-l,n + fn,n-l) ~ V^ + I(/° n+ i + /°+i,J. (114) 



Appendix C Combinatoric Derivation of Identities (|53Q and (]73p 

The ^-exponential ( |5"4T ) of the Hamiltonian kernel for the linear potential is also easy to work out 
directly, since the combinatorics in *-space are identical to the combinatorics of any associative algebra. 
In particular, the Campbell-Baker-Hausdorff expansion also holds for ^-exponentials, 

exp, (A) * exp, (B) = exp, (a + B + ±[A, B\, + 1 [A, [A, B] ^ + 1 [[A, B\, , B]^ + Cj , (115) 

Where C represents a sum of triple or more nested ^-commutators (Moyal Brackets, [A, B]+ = A*B — B-kA). 
Now, choosing A = itx and B = itp 2 + it 2 p + |it 3 , yields [A- 8 ]* = -2it 2 p - it 3 , [A, [A^IJ* = 2it 3 , 
[[A, B]^ , B}^ = 0, and hence C = 0. 

Consequently, 

exp^ (itx) * exp^ | itp 2 + it 2 p + -it 3 J = exp^ {itx + itp 2 ) . (H6) 
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But further note exp* (ax) = exp (ax), and also exp* (bp 2 + cp + d) = exp (bp 2 + cp + d). This reduces the 
^-product to a mere translation, 

exp* (ax) -k exp* [bp 2 + cp + d) = exp (ax) * exp (bp 2 + cp + d) (117) 

= exp ^ax + -iadp^j exp (6p 2 + cp + d) 

= exp ( ax + b ( p H — ia] + c ( p H — ia ] + d] 

Consequently, 

exp* (ito) * exp* ^itp 2 + it 2 p + ^it 3 ^ = exp (it (x + p 2 + i 2 /12)) , (118) 

and the identity 

exp* (it (x + p 2 )) = exp (it (x + p 2 + t 2 /l2)) (||) 

follows. 

The proof of 

exp* ( y - 2 J nhy e2X ~ Z + iVp) = ex P [~\^~'^) * e Miyp) = exp {^\ e2x ~ Z + *Vp) © 

is similar. Choosing now A = — ^e y ~ z e 2x and B = iyp, it follows that [^4, .£>]* = —2yA, so that only 
those multiple Moyal commutators survive which are linear in A. This means, then, that in the Hausdorff 
expansion [^6| for Z(A,B) = In* (exp* (^4) ★ exp* (B)), only B and terms linear in A survive. Hence, Z 
reduces to merely 

Z = B + A^ T ^_y (119) 

The Hadamard expansion in the right parenthesis means successive right ^-commutation with respect to 
B as many times as the regular power expansion of the function in the parenthesis dictates. Consequently, 



exp \--e y z e 2x j *exp(iyp) = exp* [--e y z e 2x j *exp*(%p) = exp* ( ~ ^^ e * + WPJ ■ ( 120 ) 
On the other hand, 

exp ^-^e y ~ z e 2x ^J *exp(iyp) = exp ^-I e ^+ 2 ^ exp (iy(p + i d x /2)) = exp [Af 2x ~ z + iypj , (121) 
and the identity is proven. 

Appendix D Construction of the Generating Function for the 

Liouville WFs 

From © and (|l|), it is evident that the Liouville wave functions can be generated by 

/°° dh 
I e~ lkX ^ E(k) (x), (122) 

X /sinh (tt/E^ 
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where E (k) = k 2 . Therefore, the usual wave function bilinears appearing in the WFs are generated by 
(recalling that the ij/s are real) 

exp (-e x - y cosh X) exp (-e x+y cosh Y) (123) 

l!k r <h! - lkX - iqY Mk)( x -y)M g )( x + y)- 



lsmh[7r^E{k)j °° J sinh (ny/Efa) 

Consequently, Fourier transforming this produces a generating function for WFs, 
1 

- dy e~ 2tpy exp {-e x ~ y cosh X) exp {-e x+y cosh Y) (124) 

^ J — CO 

fE(k)E(q)(x,p)- 



dk f°° dq -ikX-iqY 



lsmh[7r^E{k)J °° J sinh \ny/E~(g) 
Evaluation of this expression yields just a factor multiplying a modified Bessel function, 

/oo 
dy e~ 2ipy exp {-e x ~ y coshJf - e x+y coshY) (125) 
-oo 

™ dy e -2i P (j/+iin(coshX/coshy)) exp (_ ^ e V4 cosh X cosh coshy^ 



2 



/ r0 sh y \ ip / , \ 

I I K 2ip ( eV4 cosh X cosh y) . 



VcoshXy 

Thus, a generating function for the complete set of Liouville Wigner functions is 



2 / cosh Y\ ip / 

-( — — - | K 2 i v [ e x V4: cosh X cosh Y 

it \ cosh X J \ 



dk 



(65) 

- tkX - iqY fE(k)E( q) (x,p), 



I sinh [TTyjEik)) 00 x jsinh[TTyjE{q) 
as in the text. 

Appendix E Operator Ordering and eqn 

Given the factorized phase-space generating function 

1 „2z-2 



G(z; x,p) = exp I - - e 2x ~ z j K ip (e 2 ) , © 

what is the operator corresponding to it? According to Weyl's prescription, eq (||), the associated operator 
is 

<5(z;X,V) = T—Tn l drdadxdp Q(z; x,p) exp(ir(V — p) + ia(X — x)) (126) 
(2vr) z J 

= ^ 2 J drdadxdp exp(ir"P + iaX) exp ^— - e 2x ~ 2 — icrx^ iiQ p (e 2 ) exp(— irp). 
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The integrals over x and p may be evaluated separately, if the a contour is first shifted slightly above the 
real axis, a — > a + ie, thereby suppressing contributions to the x-integral as x — > -co. Now s = \e 2x ~ z 
gives 

j'^dx exp (-\ e 2x ~ z -i(a + ie) x) = f°° ^ (2se z )- i{a+ie)/2 exp (s) 

J~oo V 2 /JO 2s 

= 1 e -i(*+l=2)<r/2 r (_ 4 (fJ + j£ ) / 2 ) . (127) 



By (0), 



/ dpifip (e z ) exp(-irp) = i [°° dX e~ e * coshX 2^5 (X - r) = vr e~ e * coshr . (128) 



So 



0(z;^,P) = — / dTdae~ i( - z+ln2 ^/ 2 T(-i(a + ie) 12) e~ e * coshr exp(irP + iaX). (129) 

87T J 

The shifted <r contour avoids the pole in T at the origin. 

Ordering with all Vs to the right, thereby departing from Weyl ordering, yields exp(irV + ioX) = 
exp(icrA') exp (iar/2) exp(ir'P). Performing the a integration before the r integration, permits taking the 
limit e —* to obtain 

®(z;X,V) = ^Jdr (JdaY(-i((x + ie) /2) exp(iaX + iar/2 -ia(z + In 2) /2) J e~ eZcoshr exp(irP) 
= i_ ydr(47rexp(-e 2 ^-( 2+ln2 ))) e - eZcoshT exp(irP) 

= ^ f dr exp ^-Ie 2 * +r " 2 - ie 2+T - exp(irP) . (130) 



This is the operator correspondent to fl72|); it reflects the Weyl correspondence through which it was 
originally defined (although, technically, it was taken out of Weyl ordering above, merely as a matter of 
convenience, not a bona- fide change of representation). 

This form leads to a more intuitive Hilbert space representation. Acting to the right of a position 
eigen-bra, (x\ X = {x\ x, while the subsequent exponential of the momentum operator just translates, 
(x\ exp(ir'P) = (x + t\. So the full right-operation of is 

(x\<8(z;X,V) = \jdr ( x + r |exp(-ie 2 ^- 2 -ie^-i e 2 -^ 

= \ Jdy (y\ exp (-i e *+^ - \e z ^~ x - \?~^ ■ (131) 



Inserting 1 = J dx \x) (x\ gives <3(z;X,V) = f dx \x) (x\ <&(z;X,V), and leads to a coordinate space 
realization of the operator involving an x, y-symmetric kernel, 

<S(z;X,V) = ~J dxdy \x) (y\ exp ^-I e ^" 2 - le?-v+* - i e -*+i/+*j . ( 132 ) 
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This operator is diagonal on energy states: by Macdonald's identity (p7|), and the reality and orthogo- 
nality of the wave functions, 



x+y-z _ 1 z+y-x _ }_ e z-y+x 

2 2 



(E X \ <5(z;X,V) \E 2 ) = X - J dxdyr El (*) (y) exp {-^ 

= 8{E 1 -E 2 )K i ^ Wi {e z ). (133) 
This is in agreement with the corresponding phase-space expression, (|6G|). 



(134) 



The composition law of this operator also parallels its phase-space isomorph, (74), 
<8(u)<&(y) = \fdw exp (~ ( e u+v-w + e «-«+« + e - u +v+w>\ & ( w y 
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